On the nature of the optical conductivity of polaron 
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We study theoretically a fundamental issue in solids: the evolution of the optical spectra of polaron 
as the electron-phonon coupling increases. By comparing the exact results obtained by diagrammatic 
Monte Carlo method and the data got through exact diagonalization within an appropriate subspace 
of the phononic wavefunctions, the physical nature of the crossover from weak to strong coupling 
is revealed. The optical spectra are well understood by the quantum mechanical superposition of 
states with light and heavy phonon clouds corresponding to large and small polarons, respectively. 
It is also found that the strong coupling Franck-Condon regime is not accessible at realistic values 
of the coupling constant. 
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Polaron physics plays a crucial role in describing mag- 
netoresistiveperovskites high temperature supercon- 
ductors d, molecular semiconductors [1], fullerenes 
0, polar semiconductors 0], etc. However, in spite of a 
lot of efforts, there are still only a few well grounded re- 
sults concerning spectral properties, optical conductivity 
(OC) and Lehman spectral function (LSF) of polaronic 
systems at any coupling. Even the simplest model for 
periodic system, i.e. Holstein model, is not an exception 
and the nature of its spectral response at all couplings is 
still a puzzle, in particular in nearly adiabatic case, the 
most frequently encountered regime in experiments. 

In contrast to Frochlich-like models 0, which simpli- 
fies the manifold of all momenta of the lattice system 
by using continuum approximation, the Holstein model 
takes into account the true momentum space of the sys- 
tem [§]. Here the tight-binded charge carrier is coupled 
with an optical local phonon mode: 



where t is the nearest neighbor hopping amplitude, c| 
(ci) denotes electron creation (annihilation) operator in 
the site i, and &J (bi) creates (annihilates) a phonon in 
the site i with frequency ujq. The strength of electron- 
phonon interaction (EPI) is characterized by dimension- 
less coupling parameter A = g^ujo/2tD expressed in terms 
of coupling constant g and system dimensionality D. 

There are many methods able to calculate the ground 
state (GS) properties of the Holstein polaron, see Ref. Q 
as an incomplete list. In particular adiabatic and some 
variational calculations lead to discontinous transition 
between weak and strong coupling regimes as EPI in- 
creases ■ The quantum correction to this picture is a 
singular perturbation, and essentially non perturbative. 
A naive picture is that the discontinous transition be- 
comes a crossover by quantum mixture of states describ- 
ing asymptotic regimes. On the other hand, the situation 



changes drastically if some spectral properties, e.g. OC, 
are considered. The main problem is due to the infinite 
dimensional phononic Hilbert space, that has prevented 
an in-depth analysis so far. The full phonon basis exact 
diagonalization (ED) method is basically valid only 



for ID and 2D systems and suffers from effects caused by 
very small size of the considered clusters, normally lim- 
ited to 10 sites. The dynamical mean- field theory [l^ . 
though it is exact in infinite dimension system, leaves 
doubts about applicability to realistic 3D, 2D, and espe- 
cially ID systems. Finally, containing no approximation. 
Diagrammatic Monte Carlo (DMC) methods pro- 
vide approximation free results for infinite systems of any 
dimension. However, none of the above methods helps 
physical intuition or provides an understanding of the 
nature of OC and/or excited states in the weak, strong 
and especially in traditionally most puzzling intermedi- 
ate cotipling regime. An exception is momentum average 



(MA) [ij] approach. It suggests an analytic formula and 



explicitly relates the shape of OC in the strong coupling 
regime to the phononic cloud dragged by polaron mov- 
ing through the crystal. Unfortunately, the intermediate 
coupling regime, the most interesting one from an exper- 
imental point of view, is a puzzle for the MA approach 
too. 

In the present letter we present the exact DMC results 
for the OC of ID system for any interaction value rang- 
ing from the weak- to the strong-coupling regime. The 
regular part of the OC, a typical experimental probe of 
a physical system, is defined 
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Here 
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<i,j> 



Ri — R. 



(1) 



is the current operator, Ri is the position vector of the 
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FIG. 1: The empty (full) circles stand for the lattice sites 
(electron); the dashes denote phonons. The upper panel rep- 
resents a typical WCPC and the two lowest panels show pos- 
sible realizations of the SCPC. 

site i, e is charge, \ipn) and £'„ arc the exact eigenstatcs 
and eigenvalues, and a;„ = (i?„ — Eq). The lattice pa- 
rameter a and Planck constant fi are set to unity. 

To disclose the nature of the OC we introduce two con- 
trasting physical pictures of polaronic phonon cloud: i) 
weak coupling phonon cloud (WCPC), that allows a very 
good description of the weak-coupling regime; ii) strong 
coupling phonon cloud (SCPC), that is able to reproduce 
accurately the properties of the strong-coupling regime. 
The surprising result obtained by our analysis of the ex- 
act DMC data is that a straightforward linear superposi- 
tion of the weak- and strong-coupling phonon clouds, i.e. 
double phonon cloud (DPC), is capable of describing po- 
laron properties in any coupling regime. Comparison of 
the exact results with those obtained by the DPC anal- 
ysis points out that the Franck-Condon (FC) approach, 
based on the factorization of the wavefunctions associ- 
ated to charge and lattice degrees of freedom and on op- 
tical transitions in a frozen lattice, fails at all couplings 
except for strong and unrealistic values of charge-lattice 
interaction, i.e. at couplings where one can use the crud- 
est photoemission approximation (PA) that describes the 
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FIG. 2: (Color online) (a) Energy versus A in the DMC (black 
dots), DPC (dashed thick red line), lowest order perturba- 
tion [lg| (dotted magenta line), and FC (dashed thin green 
line) approaches; (b) GS spectral weight versus A in the DMC 
(black dots), DPC (dashed thick red line), WCPC (dotted ma- 
genta line), and MA(0) approaches (dashed thin green line), 
(c) Lehman spectral function in the DMC method (black dots 
for incoherent part and vertical black arrow for (5-functional 
polaron state) and in DPC approach (solid red line). 




FIG. 3: (Color online) OC in DMC method (black dots) at 
A = 0.01 (a) and A = 0.4 (b). Dashed green line is lowest 
perturbative result [3 in (a)- Dashed thin green (solid thick 
red) line in (b) is the result of WCPC (DPC) approach. 

photoionization process towards the free-electron contin- 
uum: within the FC approach one replaces the exact 
excited charge states within the GS lattice potential well 
with the free electron states in the absence of EPI. 

The phonon cloud in the weak coupling regime is not 
characterized by strong local deformations, and, thus, we 
set that the WCPC consists of up to n < JV^c phonons 
situated at arbitrary sites (upper panel in Fig. [T]). To 
the contrary, the lattice deformation is well localized in 
the strong coupling regime and, thus, we construct the 
phonon SCPC as n < Msc phonons located at 3 near- 
est neighbor sites. Due to nonlocal nature of the current 
operator Eq. ([T]) it is crucial that the deformed neighbor- 
ing sites are situated in the arbitrary place with respect 
to the charge carrier (two lower panels in Fig. [T]). We 
found that Mwc = 5 and Msc = 50 is enough to describe 
successfully the system in the range of the A values con- 
sidered in the following and concerning the ID case with 
LOo/t = 0.1 (tha lattice size is 80). 

First, we compare the GS properties obtained by DPC 
model and those coming from the exact DMC method 
and demonstrate a perfect agreement (see Fig. [5^-b). A 
very good agreement is seen even in the intermediate cou- 
pling regime suggesting that the DPC combination of the 
weak- and strong-coupling phonon clouds provides an ac- 
curate description in the whole range of coupling param- 
eters. Indeed, there are domains of A values where both 
strong- and weak-coupling phonon cloud models fail but, 
nevertheless, the results using DPC model are still in per- 
fect agreement with exact ones. Not only GS properties 
but also the LSF is well reproduced by the DPC picture 
too (Fig. One can see the obvious failure of the sim- 
plest variant of the MA approach, (i.e. MA(0) that takes 
into account only lattice deformations located at one lat- 
tice site at a time) , in describing the GS spectral weight Z 
(see Fig. [51d). This failure explicitly demonstrates that 
in the intermediate coupling regime the phonon wave- 
functions are characterized by lattice distortions involv- 
ing many sites at the same time. 
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FIG. 4: (Color online) OC in DMC (black dots), DPC (solid 
red thick line), WCPC (dashed green line), and SCPC (dotted 
magenta line) approaches at A = 0.77 (a) and A = 0.87 (b). 
Insets show comparison of OCs in SCPC and FC (solid blue 
line) approaches. The OCs in the FC approach is shifted down 
by energy « 5a;o. 



Good agreement of the OCs got by pcrturbative and 
DMC approaches at A = 0.01 (sec Fig. |3^) indicates that 
the final states of the optical transition do not contain 
more than one phonon. To the contrary, evident rise 
of OC at 2^0 for A = 0.4 (see Fig. ^jp) manifests im- 
portance of the processes where correlation between two 
successively emitted phonons is present. The agreement 
between WCPC approximation and DPC approach indi- 
cates that the contribution of the strong coupling coun- 
terpart of the phonon cloud is still negligible at A — 0.4. 

In Fig. |4]the exact data for the OC in the intermediate 
coupling regime (A « 0.8) are compared with the results 
following from the WCPC, SCPC, and DPC conceptions 
of the phonon clouds. The intermediate coupling regime 
is a stumbling stone for most of the methods because any 
approach, starting from the characteristic features of ei- 
ther weak- or strong-coupling regime, fails there. The 
results in Fig. 0] clearly reveal the reason for such diffi- 
culties, showing that the phonon cloud in this regime is 
a superposition of the WCPC and SCPC contributions. 
The dominating hump of the OC splits into two peaks: 
the former one is located at the energy determined by 
the weak coupling counterpart of the lattice deformation 
and the higher energy structure is governed by strong- 
coupling component of the phonon cloud. Indeed we 
note that if the OC obtained in the framework of the 
FC concept is shifted towards lower frequencies we re- 
cover the peak of the OC governed by strong-coupling 
lattice deformation (see insets in Fig. Hence, the OC 
in the intermediate coupling regime is not yet described 
by the FC processes where optical excitations correspond 
to electronic excitations in the rigid lattice. The high- 
est energy peak of the OC in the intermediate coupling 
regime can be interpreted as the result of light absorption 
in two successive steps: the lattice is frozen during the 
electronic transition at the first stage and then it relaxes 
and adapts itself to new electron configuration. We stud- 
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FIG. 5: (Color online) OC in DMC (black dots), DCP (solid 
red thick line), and FC (dashed blue line) approaches at A = 
1.5 (a) and A = 10 (b). Inset in (b) shows OC in FC approach 
and photoemission curve (solid orange line). Panel (b) shows 
the energy of lattice relaxation from the frozen FC state to 
relaxed excited configuration of the SCPC state. 



icd the the energy of lattice relaxation from the frozen 
FC state AE and found that it becomes negligibly small 
only in the very strong coupling regime at A > 1.2 (Fig. 
[S}:). The polaron mass at such couplings is three orders 
of magnitude larger than the bare mass. Such effective 
mass, to our best knowledge, was never observed in any 
physical system. 

By increasing A there is a net transfer of spectral 
weight towards higher frequencies, and at A = 1.5 only 
the higher energy contribution survives (see Fig. [5]a,). 
In this regime the energy of lattice relaxation from the 
frozen FC state is almost zero and the energies of the 
main peak of the OC in the DMC, DPC, SCPC, and FC 
approaches are the same. On the other hand there are 
differences in the tails of OC calculated in the above tech- 
niques, indicating the residual role played by the nonadi- 
abatic transitions even in the strong-coupling regime. Fi- 
nally, the crudest PA approximation is valid at extremely 
large couplings (see inset in Fig. [Sb)- 

Above discussed considerable success of the concept 
of the double phonon cloud paves a way to a novel all- 
coupling method called in the following Double Phonon 
Cloud ED (DPCED) approach. Limited phonon basis, 
including not the whole manifold of the phonon states but 
only the essential ones, enables studies of large systems 
up to about 100 sites. The approach is described in detail 
in the following. 

First of all we take into account the translational in- 
variance of the Hamiltonian and perform an exact diag- 
onalization of the Hamiltonian based on Lanczos algo- 
rithm, requiring that the states have a definite momen- 
tum 



17| . Each of basis vectors is a linear superposi- 



tion with appropriate phases of the translational copies 
(charge carrier and lattice configurations are together 
rigidly translated) of a state having the electron fixed 
at a site and phonon quanta located around it. 

Weak coupling counterpart is restricted to the most 5 
phonons situated up to 5 different sites which, in turn, are 
arbitrary located with respect to the charge carrier. The 
phononic Hilbert subspace corresponding to the WCPC 



4 



is generated by the following elements 
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(2) 



h=l 



Here i^jh label the lattice sites, [j] = (j/i, ft, = 1, ...5) rep- 
resents a set of 5 different sites, |0)j is the i-site phonon 
vacuum state, and the integers rij^ = 0, ....5 are such 
that "Ylih "^ih < 5. In any basis element there are up to 5 
phonons distributed at the most on 5 different sites: all 
the others sites are not deformed. In this way the scat- 
tering processes between the charge carrier and lattice up 
to 5 phonons in g-space are exactly treated. This basis 
is able to recover the self-consistent Born approximation 
and goes beyond it including vertex corrections. 

Strong coupling counterpart in the direct space repre- 
sentation contains the following set of states: 



\ph) 



(SCB) 



|0),, 



n 10). 



(3) 

Here j' indicates the nearest neighbors of j site, and nj' 
represents the number of phonons on such sites. In any 
basis element there is a cluster with at the most 3 de- 
formed nearest neighbors sites which are located at an 
arbitrary position with respect to the electron. We found 
that the sufficient basis for all studied above couplings is 
limited to {nj + J2j' n-j') ^ 50- We note that by restrict- 
ing the cluster with excited phonons to one site one comes 
to the simplest formulation of the MA approximation, i.e. 



MA(0) [18 



DPCED approach is based on the idea that the physical 
properties for any value of the EPI can be described by di- 
agonalizing the Hamiltonian in the phononic Hilbert sub- 
space generated by both above introduced sets of states 
(excluding obvious double counting). 

In conclusion, we have shown that the ground state, 
Lehman spectral function, and OC of the polaron at any 
coupling are well described in terms of superposition of 
two phonon clouds: the former one is spatially extended 
and contains small number of phonons, whereas the lat- 
ter one is well localized and includes a large number of 
vibrational quanta. Surprisingly, the lattice deforma- 
tion around polaron in the puzzling intermediate cou- 
pling regime appears to be not a peculiar characteristic of 
this particular regime but a mere superposition of quan- 
tum states representing two asymptotic couplings. The 
accuracy of the above concept of double phonon cloud 
paves the way to a novel DPCED method, which uses 
not the whole manifold of the phonon states but only the 
essential ones, and enables study of large systems which 
are inaccessible by classic exact diagonalization method. 
Comparison with the results for OC of the novel method 
with those obtained by approximation free Diagrammatic 
Monte Carlo method validates both the suggested dou- 
ble phonon cloud concept and the novel numeric method 



based on it. The novel method can be further applied 
to problems of nonequilibrium dynamics [lij where such 
methods as diagrammatic Monte Carlo fail. 
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